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ABSTRACT 

We present infrared photometry of the Galactic Bulge X-ray binary systems GX13+1 and GX17+2 
obtained in 1997 July and August using OSIRIS on the 1.8m Perkins Telescope at Lowell Observatory. 
GX13+1 clearly varies over ^0.6 magnitudes in the i^-band. Our light curve suggests a modulation on 
a timescale of ~20 days, which is in agreement with previously proposed orbital periods for the system. 
The IR counterpart of GX17+2 is also variable in the Jf-band over ~0.8 magnitudes on a timescale of 
days to weeks, extending the variability first seen by Naylor, Charles, & Longmore (1991). We discuss 
the implications our data have for Deutsch et al. 's (1999) identification of "star A" as the true IR 
counterpart of GX17+2. The variability observed in our photometry of the blend of star A and the 
foreground star NP Ser implies a ~4 magnitude intrinsic variation in the .fT-band for GX17+2. 

Subject headings: binaries: close - infrared: stars - X-rays: stars - stars: individual: GX13+1, 
GX17+2 - accretion, accretion discs 



1. INTRODUCTION 

In low-mass X-ray binaries (LMXBs), mass is trans- 
ferred from a late-type star to its highly compact com- 
panion, either a neutron star or a black hole, via an ac- 
cretion disc. These systems can be classified according 
to location within the Galaxy, accretion characteristics, 
and luminosity (van Paradijs & McClintock 1995). Map- 
ping of the X-ray sky has revealed a dozen bright X-ray 
sources within 15° longitude and 2° latitude of the Galac- 
tic Centre (Warwick et al. 1988). Known as the "galactic 
bulge" or "bright bulge" sources (henceforth referred to as 
"GBS"), these LMXBs are among the most luminous X- 
ray sources in the Galaxy (typical Lx ~ 10 38 ergs -1 ). 
However, the GBS remain the most poorly understood 
group of LMXBs, due to the heavy obscuration in the di- 
rection of the Galactic bulge which makes optical study 
nearly impossible. The GBS are thought to be neutron 
star LMXBs, and quasi-periodic oscillations (QPOs) have 
been detected in several systems (e.g. van der Klis 1989). 
Attempts to detect orbital variability in the GBS have 
generally been unsuccessful, suggesting that their periods 
may be longer than those of canonical LMXBs (Charles 
& Naylor 1992, hereafter CN92). On the basis of their 
X-ray colour-colour diagrams, neutron star LMXBs have 
also been divided into two classes, known as Z and atoll 
sources (see van der Klis 1995 for a review). In this scheme, 
six of the brightest neutron star LMXBs (which includes 
several of the GBS) are classified as Z sources, while the 
remainder fall into the atoll category. However, recently 
this distinction has come into question, as several atoll 
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sources show Z-type colour-colour diagrams when exam- 
ined on a long timescale (Muno et al. 2001, Gierlinski 
& Done 2001). Nevertheless clear differences remain be- 
tween the X-ray luminosity and spectral evolution of the 
canonical Z sources and that of the atolls. 

The infrared provides us with an ideal window for ob- 
serving these highly obscured systems. Observing in the 
IR has two primary advantages: the late-type secondaries 
in LMXBs are brighter relative to the accretion discs, and, 
more importantly for the GBS, the ratio of V- to if-band 
extinction is nearly 10 (Naylor, Charles, & Longmore 1991, 
hereafter NCL91). Over the past eight years, we have de- 
veloped a program of IR observations of X-ray binaries 
(XRBs), beginning with the discovery via colours or vari- 
ability of candidates for the IR counterparts to the X-ray 
sources using precise X-ray and radio locations (NCL91). 
Following this photometric survey we performed a spectro- 
scopic survey of LMXBs, obtaining IR spectra of a num- 
ber of systems, including Sco X-l, GX13+1, and Sco X-2 
(Bandyopadhyay et al. 1997 and 1999, hereafter B97 and 
B99). The spectra enabled us to place constraints on the 
spectral types of the mass donors in these systems, indi- 
cating that the secondaries in the GBS may be evolved 
rather than main-sequence stars. 

If, as the evidence suggests, the companion stars in the 
GBS are indeed evolved stars, we would expect any vari- 
ability in the IR light curve due to the orbital period to 
occur on a timescale of days or weeks, rather than hours 
as in canonical LMXBs. To search for long orbital periods 
in the atoll LMXB GX13+1 and the Z LMXB GX17+2 
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(as classified by Hasinger & van der Klis 1989), in 1997 
we obtained IR photometry of both sources over a period 
of approximately 6 weeks. An IR study of the GX13+1 
field mapped by the Einstein HRI in X-rays revealed a 
variable IR source at the Grindlay & Seaquist (1986) ra- 
dio position of GX13+1 (CN92). The observed variabil- 
ity made this IR source a very strong candidate for the 
counterpart to GX13+1. The i-T-band spectrum we ob- 
tained in 1997 showed Brackett 7 emission, confirming the 
identification of the IR counterpart (B97). An additional 
spectrum obtained in 1999 with improved S/N exhibited 
both Brackett 7 emission and CO absorption bands; using 
the technique of optimal subtraction, we determined the 
most likely spectral type of the secondary to be a K5 III 
(B99). 

In contrast to GX13+1, detection of the true IR coun- 
terpart to GX17+2 has proven more difficult. GX17+2 
was optically "identified" on positional grounds alone with 
a G star now known as NP Ser more than 25 years ago 
(Tarenghi & Rcina 1972). However, the absence of any 
optical variability or spectroscopic peculiarity in a star 
purported to be associated with one of the most luminous 
X-ray sources in the Galaxy has long been a major puzzle, 
with most explanations requiring that the G star be sim- 
ply a line-of-sight object. Furthermore, the IR counterpart 
to NP Ser has shown variability, and a consistent fit for 
the optical and IR colours, extinction, distance, and spec- 
tral type could not be found (NCL91). The situation has 
recently been clarified by HST observations of GX17+2 
which resolved NP Ser into two components, the fainter of 
which, labelled "star A" by Deutsch et al. (1999), has been 
suggested as the true IR counterpart to GX17+2. The un- 
usual nature of star A has been demonstrated by Callanan 
et al. 's (1999) Keck photometry which suggests that star 
A exhibits extremely large (^3.5-4 mag) IR variability. In 
this paper we present long-term (~5-week) light curves of 
GX13+1 and GX17+2 from July and August 1997; both 
sources show clear variability of greater than 0.5 magni- 
tudes in the if-band. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained if-band photometry during 14 nights in 
1997 between 16 July and 23 August using the Ohio State 
Infrared Imager/Spectrometer (OSIRIS) on the Perkins 
1.8m telescope of the Ohio State and Ohio Wesleyan 
Universities at Lowell Observatory in Flagstaff, Arizona. 
OSIRIS uses a 256x256 pixel HgCdTe NICMOS III de- 
tector; see DePoy et al. (1993) for further details on the 
design of this instrument. The images were taken with 
the f/7 camera which provides a 2.7 arcmin 2 field of view 
at a resolution of 0.63" per pixel. The summer monsoon 
season affected the observing run, causing some variabil- 
ity in the data quality; however, typical images were ob- 
tained with < 2" seeing. On each night, a series of five 
consecutive 60 second (2 second on-chip integrations x 30 
coadds) exposures were obtained for each target, with off- 
sets of ~ 20" between each exposure. The position of the 
object on the array was moved between exposures so that 
the group could later be median-stacked to produce a sky 
frame. 

The images were processed by first running a simple 
interpolation program to remove bad pixels. A median- 
combined sky image was then created from the five ex- 



posures and subtracted from each image. The resulting 
sky-subtracted frames were flat-fielded with an image con- 
structed from the difference of dome flats obtained with 
the flat-field lamps on and off. The data were then ana- 
lyzed using the IRAF reduction task "apphot" with a 4.5 
pixel aperture which was chosen to maximize the target 
signal relative to the sky background. Relative photome- 
try was performed using several field standards. Finally, 
absolute photometry of the local standards was performed. 
A journal of the IR observations with the calculated K 
magnitudes of GX13+1 and GX17+2 is presented in Ta- 
ble 1. 

3. RESULTS 
3.1. GX13+1 

Our i\-band light curve of GX13+1 is presented in Fig- 
ure 1; these data constitute the longest IR light curve of 
this source published to date. Substantial variability of 
~0.6 mag is clearly present over the six weeks of observa- 
tion. Visual inspection of the GX13+1 light curve suggests 
a characteristic variability timescale of order ~20 d. Un- 
fortunately, the IR data are not sufficiently well sampled, 
nor is the observational baseline long enough, for a formal 
statistical analysis to yield meaningful results as to the 
presence of ~20 d periodicities. The lower panel of Fig- 
ure 1 shows the ASM 1-day average X-ray light curve of 
GX13+1 during the interval of our observational program. 

The most striking feature of the CN92 light curve (which 
spanned 10 days) was an increase in brightness of ~0.4 
mag over two days (from K ~ 12.3 to 11.9). Subsequent 
nights did not show similar variability, and monitoring of 
GX13+1 over a period of 5 hrs did not show any modula- 
tion at a level greater than 0.1 mag. Groot et al. (1996; 
hereafter G96) obtained Gunn z band (0.9-1.0 yum) pho- 
tometry of GX13+1 over a period of 18 days, and found 
clear variability on a timescale of ~10 days; using a Lomb- 
Scargle periodigram they estimate the maximum power to 
be at 12.6 days, with an estimated error of ~1 day. How- 
ever, they note that the time interval covered by their 
observations was too short to formally establish a period- 
icity in the system. Analysis of 1 year of RXTE ASM 
monitoring of GX13+1 revealed a possible 25.2 d modula- 
tion in the X-ray light curve (Corbet 1996). In contrast to 
the CN92 and G96 observations, Wachter (1996; hereafter 
W96) found a 19.5 hr modulation from 4 nights of CTIO 
i-T-band photometry. They also report the detection of the 
~0.4 mag variability seen by CN92, and suggest that this 
variation occurs on a much shorter timescale than the 12.5 
or 25.2 d periods which had been suggested by G96 and 
Corbet. 

To search for evidence of our ^20 day modulation in 
the X-ray light curve of GX13+1, we re-analyzed the data 
from the RXTE ASM. The ASM light curve now spans ap- 
proximately 5 years, in contrast to the ~1 year available 
for Corbet's 1996 analysis. Using the publically available 
"definitive" ASM dwell data from the RXTE archive, we 
used the "perdgrm" task in FTOOLS v5.0.1 to perform a 
power spectral analysis on the GX13+1 data. The maxi- 
mum peak is found at ~26.7 days; however, this peak is not 
stable. When the X-ray light curve is divided into 5 ~1- 
year intervals, the frequency of the maximum peak varies 
between ~24-45 days. To illustrate the significance and 
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Table 1 

1997 K magnitudes of GX13+1 and GX17+2 



UT Date GX13+1 GX17+2 



16 July 


12 


.12±0. 


.03 


14. 


,63±0 


.25 




11 


.87±0. 


.03 


14. 


79±0 


.32 


17 July 


12 


.06±0. 


.03 








18 July 


12 


.23±0. 


.04 


14. 


,69±0 


.29 




12 


,30±0. 


.04 








19 July 


12 


,30±0. 


06 


14, 


96±0 


.38 




12 


.60±0. 


.05 


15, 


,27±0 


.43 




12 


.45±0. 


.04 








23 July 


12 


.47±0. 


.03 


14. 


,48±0 


.14 




12 


.46±0. 


.04 


14, 


,26±0 


.16 


24 July 


12 


.40±0. 


03 


14, 


,22±0 


.14 


25 July 


12 


.46±0. 


,03 


14, 


,02±0 


.12 


26 July 


12 


.30±0. 


.05 


14, 


,26±0 


.28 




12 


,24±0. 


.03 


14, 


,07±0. 


.17 


1 August 


11 


.98±0. 


.02 


14, 


30±0 


.20 


12 August 


12 


.22±0. 


.03 


14, 


83±0. 


.26 




12 


11±0. 


.03 


14, 


,81±0 


.31 




12 


,21±0. 


.03 


14, 


,87±0 


.27 


13 August 








14. 


,08±0 


.26 


14 August 


12 


.54±0. 


.04 


14, 


63±0 


.26 


20 August 


12 


10±0. 


.03 


14. 


,14±0, 


.10 


23 August 


11 


.97±0. 


.02 









variability of this modulation, we present a 2-D "dynami- 
cal" power spectrum of the GX13+1 ASM data in Figure 
2. This power spectrum is constructed from the 1-day av- 
erage light curve, summing the 3 energy channels for an 
energy range of 1.5-12 keV. In order to search for any peri- 
odic phenomena which might be transitory in nature, the 
~1700 days of data were divided into 200 day subsets, but 
with a spacing between them of only 100 days (i.e. over- 
lapping the subsets each side). Each subset was subjected 
to a Lomb-Scargle periodogram analysis covering the pe- 
riod range 4-50 days, and the resulting periodogram high- 
lights any changes in temporal behaviour. Figure 2 shows 
that variations on timescales of 20-30 days are present in 
the power spectra of GX13+1. Specifically, during the in- 
tervals of MJD 51300-51400 and 51700-51800 significant 
peaks (contour levels geq\2, corresponding to a 99% con- 
fidence level) are seen at frequencies of ~27 and ~22 days 
respectively. However, these are of a quasi-periodic nature 
and no truly periodic signal is consistently present at this 
frequency for the duration of the ASM dataset. It seems 
likely that the 25.2 d variability reported by Corbet (1996) 
from the first year of RXTE monitoring of GX13+1 is also 
part of this long-timcscale quasi-periodicity. 

On the basis of the currently available IR photome- 
try, we are not able to identify the ^20 d variability as 
an orbital modulation. A period of this length might 
also be associated with the precession of a warped accre- 
tion disc, with the true orbital period being substantially 
shorter. Cyg X-2 is an example of this scenario; it has a 
9.8 d orbital period, determined from optical spectroscopy 
(Casares, Charles, & Kuulkers 1998), and a 78 d period 
determined from X-ray light curves and attributed to disc 
precession (Wijnands, Kuulkers, & Smale 1996). If the 



ratio Pdisc/Porb is similar for GX13+1, we might expect 
an orbital period of ~2.5 days for GX13+1 for a 20 d disc 
precession period. 

However, iT-band spectroscopy has determined the most 
likely spectral type of the secondary in GX13+1 to be a 
K5iii (B99). This spectral type correlates well with a 20 d 
orbital period for GX13+1 and rules out a period of cither 
~2.5 days or the 19.5 hr period proposed by W96. It there- 
fore seems more likely that the observed ~20 d variability 
may be an orbital modulation rather than being caused by 
a precessing disc. We note that the large IR variability and 
the low Lx/L opt ratio in GX13+1 indicates that the disc 
is large and that X-ray heating of the disc and secondary 
contributes significantly to the IR flux from the system 
(CN92). This conclusion was supported by the IR spec- 
troscopy, which indicated that the K5iii secondary proba- 
bly contributes less than half of the i-C-band flux. We also 
note that the IR and X-ray light curves shown in Figure 
1 appear generally uncorrelated on a timescale of ~weeks, 
although the sampling of the IR light curve is insufficient 
for a formal correlation analysis. Therefore we suggest 
that the IR variability of GX13+1 arises from a combina- 
tion of causes, with short-term IR emission from the X-ray 
heating superimposed on an underlying long orbital modu- 
lation. This scenario would account for both the long-term 
(>10 days) variability shown in Figure 1 and the shorter 
timescale variability seen by CN92 and W96, while simul- 
taneously being consistent with the classification of the 
companion as a K5lll star (M ~ 5M Q , R ~ 25R Q ). 

3.2. GX17+2 

We observed the field of the GBS GX17+2 on 13 nights 
(see Table 1) in 1997; the if-band light curve is presented 
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in Figure 3. Variability over ^0.8 mag is clearly present 
over the six weeks of observation. Our if-band light curve 
again confirms the variability seen by NCL91, and extends 
it to show that the IR source varies over ~0.6-0.7 mag on 
a timescale of ~10 days. However, our data are not suffi- 
ciently well-sampled to determine a periodicity in the IR 
photometric modulation. To search for an X-ray period of 
comparable length, we performed a power spectral analy- 
sis of the RXTE ASM dwell data of GX17+2. Similarly to 
our analysis for GX13+1, in Figure 4 we present a dynam- 
ical power spectrum of the GX17+2 ASM 1-day average 
light curve. We find no significant peak below ^200 days. 

Using HST IR observations of the GX17+2 field, 
Deutsch et al. (1999; hereafter D99) were able to iden- 
tify a new candidate counterpart to the X-ray source, lo- 
cated 0.9" north of NP Ser. They found a photometric 
magnitude of H « 19.8 and an iJ-band reddening of 2 
for this new candidate, labelled "star A" . Subsequently, 
Callanan et al. (1999) obtained a Keck i-T-band image of 
GX17+2 and were able to resolve NP Ser and star A, de- 
riving K = 14.5 and 14.9 respectively. Archival if-band 
images from 1998 examined by Anderson et al. (1999) 
also showed star A to be only slightly fainter than NP Ser. 
However, the iJ-band magnitude and reddening found by 
D99 implies K ~18.5 for star A at the time of their obser- 
vations. Taken together, these observations suggest that 
star A exhibits ^3.5-4 magnitudes of variability in the K 
band. This extreme modulation is consistent with the pre- 
diction of D99 that the optical/IR counterpart to GX17+2 
should be highly variable. 

Our if -band images obtained with the 1.8m telescope 
at Lowell are of insufficient resolution to resolve NP Ser 
and star A. However, the blended light curve of the two 
stars, shown in Figure 3, does exhibit significant variabil- 
ity. Assuming that at the peak of its brightness, star A has 
K ~ 14.9, and that at minimum K ~ 18.5, we can esti- 
mate the effect of the variability of star A on the blended 
light curve of the two stars. At star A's minimum, we 
expect to see flux only from NP Ser at its presumably 
constant magnitude, K ssl4.5. When star A is at its max- 
imum, we observe the combined flux of both NP Ser and 
star A, which are then of roughly similar magnitude. At 
the maximum of GX17+2, then, we would expect to see 
an increase in the blended light curve of 0.75 magnitudes 
from its minimum of 14.5. This estimate is clearly con- 
sistent with the modulation shown in our GX17+2 light 
curve, which varies by ^0.8 magnitudes. Note that the 
minimum values we observed for the NP Ser/star A blend 
are consistent within the errors with K ~14.5, the mag- 
nitude of NP Ser. Therefore the photometric modulation 
seen in our blended light curve of NP Ser and star A is 



almost certainly a result of the large variability intrinsic 
to star A. The combination of the revised positional esti- 
mate for GX17+2 and the observed extreme photometric 
variability of star A are strong evidence for the identifi- 
cation of star A as the IR counterpart to GX17+2. We 
note that the ASM light curve of GX17+2 (shown in the 
lower panel of Figure 3) is remarkable for its constancy 
compared to the dramatic IR variability during this same 
interval, similar to the behaviour seen in GX13+1 (Figure 

!)■ 

4. CONCLUSIONS 

Our if-band photometry of GX13+1 shows evidence for 
a ^20 d modulation. We note that a 20 d period correlates 
extremely well with that expected for a Roche-lobe filling 
binary with a K5iii secondary (B99). On the basis of the 
observed variability in the IR together with the identifica- 
tion of the secondary's spectral type, we suggest that the 
observed ~20 d modulation may be the orbital period of 
GX13+1. We note that to conclusively identify any period 
of order ~20 d in GX13+1, it will be necessary to obtain a 
well-sampled IR light curve with a baseline of > 60 d. We 
also find evidence in the GX13+1 ASM X-ray light curve 
for quasi-periodicity on a timescale of 20-30 days, but this 
modulation is not consistently present throughout the ~5- 
year dataset. 

We also observe a ~10 d modulation in the IR light curve 
of GX17+2. Our photometry was performed on the blend 
of the foreground star NP Ser and the recently identified 
counterpart candidate star A; the observed variability is 
consistent with that expected as a result of a large modu- 
lation of the fainter of the two stars. The IR photometric 
modulation and revised position for GX17+2 strongly sup- 
ports the identification of star A as the true counterpart 
to the X-ray source. Previous attempts at obtaining an IR 
spectrum of GX17+2 have been unsuccessful, no doubt as 
a result of the small separation between NP Ser and star 
A as well as star A's large photometric variability which 
would make it nearly undetectable for a significant fraction 
of its IR period. However, now that these difficulties are 
recognized, future IR spectroscopic observations of star A 
will hopefully prove more successful in investigating the 
physical nature of the extraordinary variability of the true 
IR counterpart. 
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Fig. 1. — Upper panel: i<"-band relative photometry of GX13+1 (filled circles) from July-August 1997. A local standard is plotted for 
comparison (open triangles). Note the clear variability of GX13+1 suggesting a modulation on the order of ~20 days. Typical photometric 
errors for GX13+1 are ±0.03 mag. Lower panel: Simultaneous RXTE ASM 1-day average X-ray light curve of GX13+1. 
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Fig. 2. — Dynamical power spectrum of the GX13+1 RXTE ASM X-ray light curve. Note the significant, variable peaks around MJD 
51300-51400 and MJD 51700-51800 (contour levels >10) at frequencies of ~27 and ~22 days respectively. 
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Fig. 3. — Upper panel: K-band relative photometry of the NP Ser/star A blend, the stars at the position of GX17+2 (filled circles). A 
local standard is plotted for comparison (open triangles). Despite the large errors (generally ±0.3 mag), the blend of NP Ser and star A does 
exhibit significant variability of order ~0.8 mag. Lower panel: Simultaneous RXTE ASM 1-day average X-ray light curve of GX17+2. 
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Fig. 4. — Dynamical power spectrum of the GX17+2 RXTE ASM X-ray light curve. The highest contour level is <8, indicating that no 
highly significant peaks were found. 



